Dynamics of Large-Scale Plastic Deformation 
and the Necking Instability in Amorphous Solids 
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We use the shear transformation zone (STZ) theory of dynamic plasticity to study the necking 
instability in a two-dimensional strip of amorphous solid. Our Eulerian description of large-scale 
deformation allows us to follow the instability far into the nonlinear regime. We find a strong 
rate dependence; the higher the applied strain rate, the further the strip extends before the onset of 
instability. The material hardens outside the necking region, but the description of plastic flow within 
the neck is distinctly different from that of conventional time-independent theories of plasticity. 

PACS numbers: 62.20.Fe, 46.35.+Z, 83.60.Df, 46.05. +b 
Keywords: plasticity, necking, numerics, STZ 



-a 
o 

(N 
> 

m 
\o 
m 

o 

(N 
O 

-i— > 
c3 



i 

C 

O 

o 



X 



Conventional descriptions of plastic deformation in 
solids consist of phenomenological rules of behavior, with 
qualitative distinctions between time-independent and 
time-dependent properties, and sharply defined yield cri- 
teria. Plasticity, however, is an intrinsically dynamic 
phenomenon. Practical theories of plasticity should con- 
sist - not of intricate sets of rules - but of equations 
of motion for material velocities, stress fields, and other 
variables that might characterize internal states of solids. 
Roughly speaking, a theory of plasticity, especially for 
an amorphous solid, should resemble the Navier-Stokes 
equation for a fluid, with the pressure replaced by a stress 
tensor, and the viscous forces replaced by a constitu- 
tive law relating the rate of plastic deformation to the 
stresses and internal state variables. That constitutive 
law should contain phenomenological constants, analo- 
gous to the bulk and shear viscosities, that are measur- 
able and, in principle, computable from molecular the- 
ories. Yield criteria, work hardening, hysteretic effects, 
and the like would emerge naturally in such a formula- 
tion. 

The goal of the STZ (shear-transformation-zone) the- 
ory of plasticity |J ||, ||, |], ||, [| 0, from its inception, 
has been to carry out the above program. In this paper 
we show how the STZ theory describes a special case of 
large-scale yielding, specifically, the necking instability of 
a strip of material subject to tensile loading. There is a 
large literature on the necking problem. References that 
we have found particularly valuable include papers by 
Hutchinson and Neale ||, McMeeking and Rice and 
Tvergaard and Needleman ]To| . Our purpose here is to 
explore possibilities for using the STZ theory to investi- 
gate a range of failure mechanisms in amorphous solids, 
possibly including fracture. We are able to follow the 
necking instability far into the nonlinear regime where 
the neck appears to be approaching plastic failure while 
the outer regions of the strip become hardened and re- 



main intact. We find that necking in the STZ theory is 
rate dependent; the instability occurs at smaller strains 
when the strip is loaded slowly. One especially important 
element of our analysis is our ability to interpret flow and 
hardening in terms of the internal STZ variables. 

To make this problem as simple as possible, we consider 
here only strictly two-dimensional, amorphous materials. 
By "strictly," we mean that elastic and plastic displace- 
ment rates are separately planar as in two-dimensional 
molecular dynamics simulations. The two-dimensional 
STZ-equations presented in this paper are based on ear- 
lier work by Falk, Langer, and Pechenik ^ |ll|, We 
use Eulerian coordinates in which, as in fluid dynamics, 
the variables Xi denote the current physical positions of 
material elements. Let the system lie in the x\ = x, 
x 2 = y plane, and write the stress tensor in the form: 



= -p b~i. 



p = — \ (Jkk , where p is the pressure 



and Sij is the deviatoric stress - a traceless, symmetric 
tensor. In analogy to fluid dynamics, let Vi{x,y,t) de- 
note the material velocity at the physical position x, y 
and time t. Then the acceleration equation is: |13|] 
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Here, p is the density which, because we shall assume a 
very small elastic compressibility and volume conserving 
plasticity, we shall take to be a constant. The symbol 
d/dt denotes the material time derivative acting on a 
scalar or a vector field: 



d _ d 
~dt = di 
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Our first main assumption is that the rate of deforma- 
tion tensor can be written as the sum of linear elastic and 
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plastic contributions: 
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where [i is the shear modulus, K = + i>*)/(\ — v*) 
is the two-dimensional inverse compressibility (or bulk 
modulus), and v* is the two-dimensional Poisson ratio. 
The symbol T>/T>t denotes the material time derivative 
acting on any tensor, say Aif. 
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and 10 ij is the spin: 
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The plastic part of the rate-of-deformation Dfj aat , like 
s.^, is a traceless symmetric tensor, thus the plastic de- 
formations are area-conserving. For present purposes, 
we use a simple, quasi-linear form of the STZ theory in 
which 



Df- ast = e qij(s, A); 



Qij(s, A) 



(6) 



and eo is a material-specific constant. The traceless, sym- 
metric tensor Aij is the internal state variable mentioned 
earlier. It is proportional to a director matrix that spec- 
ifies the orientation of the STZs; its magnitude is a mea- 
sure of the degree of their alignment. The equation of 
motion for is: 



VAjj _ _ 1 



A,, 



(7) 



In Eq. (|6|), A plays - very roughly - the role of the "back 
stress" or "hardening" parameter in conventional theories 
of plasticity 0, |l5|, [l6|, a major difference being that A 
emerges directly from a rate equation governing the pop- 
ulation of STZs and is, in principle, a directly measurable 
quantity jl], [l7j . If the second term on the right-hand side 
of Eq. (|7) were missing, then A would be proportional 
to the integrated plastic strain. This second term, how- 
ever, which is produced by the creation and annihilation 
of STZs, is a crucial element of the STZ theory. As we 
shall show briefly below, this term produces the exchange 
of dynamic stability between viscoelastic and viscoplas- 
tic states that replaces the conventional assumptions of 
yield surfaces and other purely phenomenological rules of 
behavior. 

With one important exception, Eqs. (||) and (Q) con- 
stitute a tensorial version of the original STZ theory 
obtained by linearizing the stress dependence of the 
rate factors and rescaling. Because of the linearization, 



these equations do not properly describe memory effects 
present in the full theory that are important when the 
system is unloaded or reloaded, but this will not affect 
our results until the system reaches the necking instabil- 
ity. Only after the neck starts to flow plastically, causing 
the hardened regions to unload, will we need the full 
non-linear theory to determine if the observed behavior 
is pertinent. We have chosen the rescaling so that all 
stresses and moduli are expressed in units of the plastic 
yield stress. We also have assumed that the local density 
of STZs is always at its equilibrium value so that we do 
not need to solve an extra equation of motion for that 
field (denoted by the symbol A in earlier papers). 

The important exception alluded to above is the pres- 
ence of the absolute- value bars in Eq. (|^) . The expression 
inside the bars is proportional to the rate at which plastic 
work is being done on the system, a quantity which ap- 
pears in the original theory as a non-negative factor in the 
STZ annihilation and creation rates. A negative value of 
this quantity would be unphysical. In earlier studies of 
spatially uniform systems, this quantity always remained 
positive; however, we have observed negative values in 
the present calculations. The absolute value prevents 
such unphysical behavior and is consistent with the in- 
tent of the original theory. We emphasize, however, that 
this term contains some of the principal assumptions of 
the STZ theory. There are other possibilities for it (see, 
for example, S) and, as yet, there is no first-principles 
derivation. 

To understand the transition between viscoelastic and 
viscoplastic behaviors at the yield stress, and the role 
played by the state variable A, it is easiest to look first at 
a uniform system under pure shear. Let s xx = — s yy = s, 
Sxy = 0, A xx = —A yy = A, A xy = 0] and consider a 
situation in which s is held constant. Eqs. (||) and (Q) 
become 



e = e-o (s - A) 
A = (s-A)(l-sA), 



(8) 
(9) 



where e is the total strain rate. At s = 1 , these equations 
exhibit an exchange of stability between the non-flowing 
steady-state solution with e = 0, A = s for s < 1 and 
the flowing solution with e ^ Q, A = 1/s for s > 1. As 
explained in earlier publications, the steady-state system 
is "jammed" or "hardened" in the direction of the applied 
stress for s < 1 ; whereas, for s > 1, new STZs are being 
created as fast as existing ones transform, and there is a 
nonzero plastic strain rate. 

Our goal now is to see how this exchange of stabil- 
ity occurs in a dynamic, spatially nonuniform situation. 
Consider a rectangle with straight grips at x = ±L(t). 
The upper and lower surfaces, at y — ±Y(x,t), are free 
boundaries. We assume symmetry about both the x and 
y axes so that we need to consider only the first quadrant 
of the system. On the free upper boundary, the relation 
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between the material velocities and the motion of the 
surface is 



dY 



dt 



— = v y (x,Y,t)-v x (x,Y,t)[ — \ 
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We also must specify stress conditions on this surface: 
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Here, 7 is the surface tension, k is the curvature, and the 
subscripts n and t denote normal and tangential compo- 
nents respectively. The grips at x = ±L(t) move out- 
ward at a predetermined strain rate, L/L = fi; thus 
v x (L,y,t) = LVt for < y < Y(L,t). Note that we 
do not constrain v y along this edge; we allow the grip to 
slide in the y direction. 

We wish to study how the shape of the upper surface, 
Y(x, t), changes as the grips on the sides are moved out- 
ward at various strain rates f2. Rather than trying to 
track this surface through the most general possible de- 
formations, we assume that Y(x, t) remains single valued 
and simply make a change of variables: 
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We then transform Eqs. dl|), (||), and (0), and the bound- 
ary conditions (|l^) and ( |ll"|) into equations of motion for 
the velocity, the stress, the state variable A, and the mov- 
ing boundary Y , all expressed as functions of £, C an d t. 
We solve these equations in the fixed square < £ < 1, 
< C < 1. 

In all of the calculations described here, we have used 
p = 1, (J, = 100, K = 300, and 7 = 0.1. Our initial 
conditions are L(0) — 4 and Y(x,0) = 1 — 6(x), where 
S(x) = 0.01 exp(— 8a; 2 ) is a small deformation that breaks 
translational symmetry. We chose two values for e : 0-1 
(hard) and 0.3 (soft), and two values for the strain rate f2: 
0.01 (fast) and 0.001 (slow). The time taken by sound 
waves to cross the system is approximately L \f~pfp ~ 
0.4. This is smaller than the characteristic time scale for 
plastic deformation, which we have scaled to unity, and is 
much smaller than the actual time scales that we observe 
for our relatively small pulling rates f2. Thus, our system 
is elastically quasi-stationary, and the precise value of p 
is not important. 

We have solved these equations on a fixed, non-uniform 
80 x 20 grid in £, ( space, using the implicit differential- 
algebraic solver DASPK |H . In order to suppress numer- 
ical instabilities, we have added a small viscosity pi] V 2 Vi 
to the right-hand side of the acceleration equation ([l]), 
and have set 77 = 0.1. 

Fig. [l] shows initial and final shapes of samples under- 
going tensile tests for four different combinations of the 
two parameters and £q as indicated. For clarity, we 



a 
b 



FIG. 1: Initial and final shapes of the material in four nu- 
merical tensile tests: (a) eo = 0.1, fl = 0.001 (b) eo = 0.1, 
O = 0.01 (c) eo = 0.3, Q = 0.001 (d) e = 0.3, Q = 0.01. 
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a: e = 0.1, n = 0.001 

b: e = 0.1, n = 0.01 

c: e = 0.3, a = 0.001 

- d: e„ = 0.3, n = 0.01 
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FIG. 2: The engineering stress a xx at the grip plotted against 
the engineering strain i xx for the four cases shown in Fig. 
The big circle marks the state whose internal properties are 
shown in Fig. ^. 



show the complete strip although only the the behav- 
ior of the upper right quarter was computed. The final 
shapes were arbitrarily chosen at the time when the en- 
gineering stresses at the grips were roughly half of their 
peak values (see Fig. |^). There is a necking instability in 
all four cases, but it occurs at greater strain for the faster 
pulls. This rate dependence is also apparent in Fig. |^, 
which shows the engineering stress at the center of the 
grip, <r xx (L, 0, t) = a xx (L,0,t)Y(L,t)/Y(L,0) as a func- 
tion of the engineering strain e xx = [L(t) — L(0)]/L(0) for 
all four cases (remember that a xx = s — p, and s = s xx ). 
Thus, although the "softness" parameter eo controls the 
overall plastic response of the material, the onset of the 
necking instability is controlled by the applied strain rate. 

To see what is happening internally, we show in Fig. ^ 
graphs of s and A(= A^) along the centerline of the strip 
(the x-axis) for case (a) in Fig. || shortly after the sample 
is starting to neck. According to Eq. @, the plastic flow 
rate is proportional to s — A. Outside the necking region, 
s = A < 1; thus the system in this region has hardened 
and deforms only elastically. Inside the necking region, 
however, s rises well above unity and A becomes small. 
Here the system has come close to steady-state flow on 
the A = 1/s branch of stationary solutions of Eq. (^). 
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FIG. 3: The deviatoric stress s, and A, along the a;-axis. This 
picture corresponds to case (a) indicated by the big circle in 
Fig. ^ shortly after the sample is starting to neck. 

This internal structure of the STZ picture of necking 
dynamics makes it clear that the strain-rate dependence 
shown in Figs. [I] and || is caused by the competition be- 
tween the rate of elastic loading and the rate at which 
hardening occurs, the latter being governed by the equa- 
tion of motion for A, Eq. (0). When the loading is slow, 
A grows along with the stress s, and there is little plastic 
flow anywhere until the stress exceeds the yield stress in 
the necking region. In the opposite limit, when the load- 
ing is fast, A remains appreciably smaller than s for a 
longer time during which the material undergoes plastic 
deformation everywhere. It would be useful to test this 
prediction of the STZ theory by measuring necking, say, 
in amorphous metals. We presume that various ingredi- 
ents of the full STZ theory, such as stress-dependent rate 
factors and other features that have been ignored here, 
would be needed to fit experimental data quantitatively, 
and that we would learn much about the theory from 
such an effort. 

Note that the behavior shown in Fig. ||is quite different 
from that predicted by conventional, time-independent 
plasticity theory, in which there would be a plastic zone 
with s = 1 inside the neck, that is, s would remain at the 
yield stress. We have found no evidence that this con- 
ventional behavior occurs in the simulations presented 
here, even for the smallest pulling speeds. (For cavita- 
tion, the STZ theory predicts a conventional plastic zone 
around a growing hole when the growth rate is very slow 
Once the instability sets in, the development of the 
neck is governed by the elastic energy already stored in 
the strip. We have confirmed this feature of late-stage 
necking dynamics by performing numerical experiments 
in which we stop the motion of the grips, that is, hold 
them fixed, at various times after the neck has started to 
form but well before it has grown appreciably. We find 
that so long as the stored elastic energy is large enough, 
stopping the remote loading in this way has almost no 
effect on the neck; it continues to grow just as before, 
driven by the elastic unloading. 



The behavior described in the last sentence - necking 
driven by stored elastic energy - looks in many ways like 
fracture, although necking differs from ordinary fracture 
in that the stress concentration that triggers the insta- 
bility is due to narrowing of the strip as a whole rather 
than to a localized defect on just one surface. Neverthe- 
less, the behaviors shown in Fig. |^ (and other, later-stage 
results not shown here) suggest the onset of a localized, 
propagating failure mechanism. In order to study the 
connection between necking and fracture in adequate de- 
tail we believe that we shall need to use the full STZ 
theory and to improve our numerical resolution. 
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